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Abstract. The field of Diffraction Dissociation, which is the subject of this workshop, began 50 
years ago with the analysis of deuteron stripping in low energy collisions with nuclei. We return to 
the subject in a modern context- deuteron dissociation in ^/snn — 200 Ge V d-Au collisions recorded 
during the 2003 RHIC run in the PHENIX experiment. At RHIC energy, d^n+p proceeds predom- 
inantly (90%) through Electromagnetic Dissociation and the remaining fraction via the hadronic 
shadowing described by Glauber. Since the dissociation cross section has a small theoretical error 
we adopt this process to normalize other cross sections measured in RHIC. 



^ : INTRODUCTION 

<N . 

When deuteron beams were first accelerated to 190 MeV and collided with internal tar- 
gets in the Berkeley cyclotron, experiments found a very collimated forward beam of 
neutrons which were identified with the process of absorptive stripping originally pro- 
posed by Serber[l]. Glauber[2] then showed that deuteron breakup could also proceed 
via a process he called "free dissociation", which has no classical analog. Absorption of 
part of the deuteron occurs even when neither nucleon strikes the target nucleus (treated 
as a black disc) and this absorption mixes unbound states of the proton and neutron 
which can then appear in the outgoing beam. The calculated cross section for this pro- 
cess is large (60% of the absorptive stripping cross section). 

At RHIC collision energy (^snn = 200 GeV) a second mechanism for free dissocia- 
^ ■ tion of the deuteron becomes dominant. The intense Coulomb field of the target nucleus 

appears to the incident deuteron as a beam of photons whose flux can be calculated 
by the Equivalent Photon Approximation originally due to Fermi[9, 8]. Since the spec- 
trum extends well above the deuteron photodissociation energy (Ey = 2.23 MeV) this 
becomes the dominant process for free dissociation. A recent calculation[5] for RHIC 
yields 1.38 (±5%) barn of which 0.14 barn is due to the original nuclear dissociation 
process calculated by Glauber. 

In 2 companion papers in the early 60's Good and Walker[7] observed that both 
Coulomb dissociation and nuclear dissociation should have analogs in diffractive excita- 
tion of elementary particles in high energy collisions with nuclei. Diffractive processes 
play a significant role in both Coulomb interactions (in Ultraperipheral Collisions) and 
hadronic collisions of Heavy Ions. This should be contrasted with e-p (HERA) and pp 
(Tevatron) where either one or the other process is studied. A unique aspect of the Heavy 
Ion program in ATLAS is that HERA and Tevatron measurements can be extended to 
higher energies and nuclear targets [10, 11]. 
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D-AU CROSS SECTIONS 

In addition to the dissociation cross section the total d-Au inelastic cross section is of 
interest for the RHIC program. The inelastic cross section is sampled in the experiments 
("min-bias trigger") for use as a luminosity monitor. Once the min-bias cross section 
is known those of other processes recorded during the same luminosity interval can be 
calculated in the usual way. There are 2 approaches to this cross section normalization. 
In the first it is derived from known, elementary, NN inelastic cross sections using a 
Glauber model with a Woods-Saxon distribution parametrization of the p,n distributions 
in the nuclei. Calculations done for d-Au at RHIC energy range from 2.26±0.1 barn 
[12] to 1.93 barn[6]. It is difficult to assign an overall error to this calculation since one 
may find discussions in the literature of whether or not to include the diffractive part of 
NN cross sections and whether n and p should have the same matter distributions in the 
Au nucleus, etc. 

The second approach, which we adopt here, is to directly determine the min-bias 
trigger cross section by comparison to the reliably calculated [5] deuteron dissociation 
process which was also measured in the PHENIX 2003 data. 



INSTRUMENTATION 

The four RHIC experiments have mid-rapidity spectrometers with different character- 
istics but all share identical Zero Degree Calorimeters (ZDC's) located at ±18 m. The 
ZDC's cover ±5 cm (in x and y) about the forward beam direction and have an energy 
resolution of Oe/E < 21% for 100 GeV neutrons within x,y< 4.5 cm[3, 4]. Essentially, 
all non-interaction ("spectator") neutrons are detected in the ZDC's, while charged par- 
ticles are mostly swept out of the ZDC region by strong ( 16T ■ m) accelerator dipoles at 
z=±llm. 

The same dipoles sweep spectator protons from d-dissociation beyond the outgo- 
ing beam trajectory (since they have twice the deuteron charge-to-mass ratio) and in 
PHENIX they are detected in a proton calorimeter("fCal")[14]. 

The PHENIX experiment used two additional hodoscopes (BBC's)[4], located at 
z= ±1.5m and covering 3.0 < \t\bbc\ < 3.9, as its main min-bias trigger. Events with one 
or more charged particles hitting both the +z and -z BBC fired this trigger. Determining 
Obbc is equivalent to determining the luminosity of the PHENIX data. d-Au events 
occuring well within the z= ±1.5 m interval between the BBC's fire this trigger with 
an efficiency of 88 ±4% [13] but the efficiency falls off for \z eV ent\ > 40 cm. For this 
reason we will determine Obbc using only events within this interval. A correction is then 
applied for the fraction of all RHIC events within this interval. The actual distribution 
of events along z within the data can be measured using time-of-flight measurements 
between the ZDC's for events with a ZDC coincidence (with single event resolution of 
Oz ~ 2 cm). 
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FIGURE 1. Energy deposition in proton calorimeter vs. ZDC (neutron calorimeter) for events with 
some activity in the ZDC in the deuteron beam direction. This sample includes absorptive stripping as 
well as d^ n+p. 



DATA ANALYSIS 

Typical event rates were several kHz (J??~ 1 — 4 x 10 28 cm~ 2 s~ l ) for all processes 
considered here. Therefore this analysis is based on a representative data sample with 

N B % = 230k (1) 

N trig [ZDC Au "of 'ZDC ' d ) = 460k (2) 

events, where subscripts Au and d represent Ezdc > 10 GeV in the Au or d direction. 
The second trigger is sensitive to d-dissociation, characterized by a 100 GeV neutron in 
ZDCd and a 100 GeV proton in fCal with no activity at mid-rapidity. 

Additional data samples were recorded with one of the RHIC beams intentionally 
displaced by up to 1 mm to measure the fraction of triggers due to d-Au collisions (as 
opposed to beam gas background). The largest background was <3% and quoted rates 
have been corrected for the measured background. 

The BBC rate, corrected for accelerator interaction distribution is 

N c B °Bc Cted = 228634events (±0.5%) (3) 
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FIGURE 2. Neutron impact distribution for both absorptive stripping and dissociation events. 



D-DISSOCIATION ANALYSIS 

As stated above, d-dissociation events have a clear signature in the PHENIX experiment. 
This is illustrated in Figs 1 where we display the forward (in the deuteron direction) 
neutron vs. proton calorimeter energy for ZDC-trigger events. 

The neutron impact parameter distribution (measured by the PHENIX Shower max- 
imum detector) is displayed in Fig. 2 both for events with the d-dissociation signature 
("Coulomb") and BBC trigger ("hadronic"). The latter correspond exactly to Serber's 
absorptive stripping proccess and, as noted in ref.[6], the neutron angular distribution 
should have an interesting correlation with event centrality. This will be discussed in a 
future note. 

In any case one can see that neutrons have a small angular divergence and conse- 
quently there is only a small correction for ZDC acceptance. Instead, the dominant cor- 
rection is for absorptive stripping events which feed into the dissociation sample. In 
order to extract the dissociation event yield we used an iterative procedure, fitting the 
sum of ZDC^+fCal total energy to the sum of 100 +200 GeV lineshapes and correcting 
for calculated efficiency as successive cuts on activity in other detectors were applied. 
The first 2 iterations yield 

N(d^n + p) = 157149and 156951 (4) 

events, so the procedure is clearly stable. 
Our final result is : 

a B BC = N(BBC)/N(d^n+p) x a (d^n + p) =228634/158761 x 1.38(±0.5%) (5) 
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= 1.99 (±1.6% ±5.0%) barn. 



(6) 



This is the quantity needed for luminosity normalization. 

In order to compare with Glauber calculations in the literature we then correct for the 
BBC detector efficiency given above: 

c {inelastic d - Au ) = o corrected (BBC) = 1.99/0.88 (7) 

= 2.26(±1.6%±5.0%±4.5%)barn (8) 

where the last 2 errors reflect the theoretical error on o(d — > n ± p) and the BBC 
inefficiency uncertainty. 

A similar analysis yields the cross section for the selection ZDCa u (E > IQGeV), also 
used as a min-bias trigger: 

o{ZDC Au ) = 2.06(±1.7%± 5.0%) barn. (9) 
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